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Preliminary  experimental  results  of  forced  convection  condensation  of  steam  inside  a  horizontal  microtube  are  reported.  The 
stainless  steel  microtube  has  a  diameter  of  250  |tm  and  a  length  of  48  mm.  The  inlet  saturated  vapor  Reynolds  number  ranges  from 
390  to  4  JOO  and  the  average  inlet  saturated  vapor  velocities  range  from  33  to  276  m/s.  The  mean  subcooled  temperature  varies  from 
12  to  25°C.  The  two-phase  pressure  drop  ranges  from  near  zero  to  65  kPa.  The  vapor  quality  at  the  inlet  and  outlet  are  100%  and  24 
to  94%,  respectively.  The  experimental  average  Nusselt  numbers  are  compared  to  those  from  stratified  condensation  flow  and  annular 
condensation  flow  in  macrotubes  found  in  the  literature.  The  experimental  average  Nusselt  numbers  are  (1)  3  to  50  times  smaller 
than  the  macrotube  values,  (2)  do  not  depend  on  the  inlet  saturated  vapor  Reynolds  number  or  velocity,  and  (3)  increase  nearly 
linearly  with  the  increase  in  the  change  of  vapor  quality  from  the  inlet  to  the  outlet.  The  significant  deviation  from  macroscale 
theory  observed  may  be  due  to  different  condensation  mechanisms  caused  by  the  relatively  large  droplet  sizes  relative  to  the  inside 
diameter  of  the  microtube.  The  two-phase  pressure  drop  should  include  the  contribution  of  surface  tension. 
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ABSTRACT 


This  thesis  reports  preliminary  experimental  results  of  forced  convection  con¬ 
densation  of  steam  inside  a  horizontal  microtube.  The  stainless  steel  microtube  has 
a  diameter  of  250/zm  and  a  length  of  48mm.  The  inlet  saturated  vapor  Reynolds 
numbers  range  from  390  to  4, 700,  and  the  average  inlet  saturated  vapor  velocities 
are  from  33  to  218ml s.  Dimensionless  vapor  velocities  (ratio  of  vapor  shear  axial 
gradient  to  gravity  radial  gradient)  are  from  15  to  150,  and  liquid  volumetric 
fraction  in  the  vapor  is  from  2.3  x  10”“*  to  41  x  10““*.  The  mean  subcooled 
temperature  varies  from  12  to  25"C.  Two  phase  pressure  drops  range  from  near 
zero  to  8hkPa.  The  vapor  quality  at  the  inlet  and  outlet  is  100%  and  24  ~  94%, 
respectively.  The  experimental  average  Nusselt  numbers  are  compared  to  those 
from  the  stratified  condensation  flow  and  annular  condensation  flow  in  macrotubes 
found  in  the  literature.  It  is  found  that  the  experimental  average  Nusselt  numbers 
(1)  are  3  to  50  times  smaller  than  the  macrotube  results,  (2)  do  not  depend  on  the 
inlet  saturated  vapor  Reynolds  numbers  or  inlet  saturated  vapor  velocities,  and  (3) 
increase  nearly  linearly  with  the  increase  in  the  change  of  vapor  quality  from  the 
inlet  to  outlet.  The  significant  deviations  from  macroscale  theory  observed  may  be 
due  to  a  different  condensation  mechanism  caused  by  relatively  large  droplet  sizes 
compared  to  the  inside  diameter  of  the  microtube.  The  two-phase  pressure  drop 
should  include  the  contribution  of  surface  tension. 
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CHAPTER  1 


INTRODUCTION 

In  the  past  decade  there  has  developed  a  growing  interest  in  microscale  systems 
and  devices,  or  microelectromechanical  systems  (MEMS)  technology.  Potential 
applications  of  MEMS  technology  may  be  found  in  fields  such  as  bioengineering 
and  biotechnology,  aerospace,  mini  heaters  and  electronic  cooling  heatsinks  and 
heat  exchangers,  electronics  and  microelectronics,  material  processing,  thin  film 
deposition  technology,  micro  optics,  micro  sensors  and  processors,  micro  heat  pumps 
and  heat  pipes,  and  micro  vapor  compression  refrigeration  systems.  Pa^t  work  on 
MEMS  technology  has  primarily  been  directed  toward  fabrication  techniques  and 
considerable  progress  has  been  made  in  this  area.  As  a  consequence,  emphasis  in 
MEMS  technology  is  beginning  to  include  the  analysis  and  design  of  micro  systems. 
Better  understanding  of  microscale  flow  and  heat  transfer  will  help  to  assess  the 
combined  influences  of  mechanical,  electromagnetic,  thermal,  optical,  and  fluid 
properties  and  processes  [1,  2,  3], 

Two  general  areas  of  research  concerning  microscale  flow  and  heat  transfer  are 
being  studied.  One  is  single-phase  flow  and  heat  transfer,  and  the  second  concerns 
two-phase  flow  and  heat  transfer  with/without  phase  change.  To  date,  however, 
work  in  investigating  two-phase  condensation  in  microchannels  or  microtubes  has 
not  been  found.  If  the  main  task  in  the  previous  study  of  microscale  two-phase 
flow  boiling  and  heat  transfer  has  been  directed  toward  removing  heat  from  high 
heat  flux  components  (microelectronic  chips,  LSICs,  laser  devices,  etc.),  it  was 
understandable  why  microscale  two-phase  condensation  has  not  yet  been  dealt  with. 
The  importance  of  understanding  the  mechanisms  and  fundamental  differences 
involved  with  flow  and  heat  transfer  in  two  phase  condensation  in  microgeometries 
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can  be  seen  in  some  promising  applications  including  the  design  of  industrial  mini 
condensers,  the  construction  of  micro  vapor  power  systems  and  micro  vapor  com¬ 
pression  refrigeration  systems,  and  the  design  of  micro  heat  pipes,  etc.  Therefore, 
it  is  the  lack  of  the  work  in  this  area  and  the  potential  microscale  applications  that 
provoke  the  interest  of  the  current  research. 

1.1  Review  of  Microsccde  Flow 
and  Heat  Transfer 

General  strategies  used  for  MEMS  technology  are  to  microfabricate  individual 
components  of  a  large-scale  system  in  order  to  improve  the  system,  and/or  to 
design  and  build  a  completely  new  microscale  device  with  functionality  comparable 
to  a  macroscale  system.  However,  the  downsizing  of  components  or  systems  may 
result  in  some  special  effects.  The  scaling  effects  are  one  of  the  most  important 
factors  to  be  assessed  in  the  design  process  [4].  For  example,  the  linear  increase 
in  surface-area-to-volume  ratio  with  scale  reduction  of  an  energy  conversion  and 
transport  system  will  improve  surface  processes  such  as  convection  heat  transfer 
and  adsorption  [4].  The  scaling  effect  in  convection  heat  transfer  can  be  shown  by 
an  example  of  internal  convection.  In  laminar  fully-developed  flow  (most  microtube 
flow  is  laminar)  with  characteristic  dimension  d,  the  heat  flux  scales  as  d~^  while 
the  pressure  drop  scales  as  d~*  for  constant  mass  flow  rate.  Thus,  the  heat  transfer 
is  much  less  sensitive  to  scale  than  the  pressure  drop.  Scaling  effects  for  external 
convection,  radiation,  and  phase  change  processes  are  more  complicated. 

Sizes  of  MEMS  components  or  systems  range  from  micrometers  to  millimeters. 
At  these  scales,  in  addition  to  the  scaling  effects  mentioned  above,  actual  processes 
associated  with  fluid  flow  and  heat  transfer  in  microgeometries  may  deal  with  a 
number  of  other  microscale  factors  such  as  viscosity  variation  near  the  surfaces, 
slip  flow  at  the  boundaries,  compressibility,  and  microcontinuum  effects  [5].  For 
example,  it  has  been  shown  that  for  gas  flow  in  very  small  channels  the  continuum 
assumptions  may  break  down  when  the  molecular  mean  free  path  is  on  the  order 
of  the  dimensions  of  the  channel  [6,  7,  8].  Other  observed  microscale  effects  in¬ 
cluded  two-  and  three-  dimensional  transport  eflTects  [9],  and  effects  of  temperature 
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variations  of  the  transport  fluid  [10]. 

By  far,  the  majority  of  the  reported  research  in  microscale  flow  and  heat  transfer 
is  based  on  empirical  or  semiempirical  correlations  of  experimental  results  and  no 
universally  accepted  theoretical  models  are  available. 

Single-phase  flow  and  heat  transfer  in  microgeometries  has  been  the  most  ex¬ 
tensively  investigated  field.  This  research  was  initiated  by  the  work  of  Tuckerman 
[11]  who  used  microchannels  of  about  50/xm  to  cool  micro  electronic  chips.  A 
IQOW/cm^  heat  flux  was  removed  from  the  silicon  chips  using  water  with  a  pressure 
drop  of  about  2  atm  and  a  7l°C  temperature  rise  above  the  fluid.  Since  then,  a 
considerable  number  of  related  research  projects  have  been  conducted.  Phillips  [12] 
developed  a  complete  resistance  model  to  predict  heat  transfer  and  flow  performance 
of  microchannels  by  improving  on  the  Tuckerman  model  [11].  Weisberg  et  al.  [13] 
conducted  a  numerical  analysis  of  heat  transfer  in  a  microchannel  and  acquired  the 
profile  of  the  Nusselt  number  along  the  walls  of  the  microchannel.  Yu  et  al.  [14] 
developed  a  simple  resistance  model  to  design  an  air-cooled  microchannel  heatsink 
used  for  the  cooling  of  a  high  power  Traveling  Wave  Tube  (TWT).  In  addition  to 
theoretical  investigations,  several  small  scale  devices  or  systems  that  make  use 
of  single-phase  convection  have  been  developed.  Examples  include  small  scale 
silicon  thermal  flow  sensors  [15],  microscale  pressure  sensors  [16],  micromembane 
pumps  [17],  cooling  devices  for  carbon  dioxide  laser  arrays  [18],  a  mini  heatsink 
for  a  high  power  TWT  [14],  etc.  A  maximum  heat  flux  of  IZOOW/cm^  for  water 
cooling  has  been  reported  by  Koh  et  al.  [19].  Recently,  Yu  at  al.  [14]  reported 
an  air-cooled  microchannel  mini  heatsink  that  was  used  to  remove  a  heat  flux  of 
up  to  14W/cm^  from  a  high  power  TWT  with  the  characteristic  of  low  pressure 
drop.  However,  the  models  used  for  the  above  investigations  are  mostly  based  on 
macroscale  assumptions.  For  instance,  the  Nusselt  number  in  the  microchannel 
usually  takes  on  a  value  that  is  from  that  acquired  between  two  macroscale  infinite 
parallel  flat  plates  with  an  isoflux  boundary  condition  when  the  aspect  ratio  of  the 
microchannel  is  large  enough  [19,  12,  14,  20].  Such  a  treatment  would  undoubtedly 
ofiFer  a  means  to  realize  preliminary  engineering  applications,  but  it  may  have 
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obscured  possible  microscale  effects  due  to  the  approximations  used  in  the  models 
[14,  20].  Although  a  large  number  of  studies  have  been  directed  toward  detecting 
microscale  effects  for  single  phase  flow  and  heat  transfer,  the  current  research 
seems  not  to  be  able  to  reach  conclusive  results  due  to  the  insufficiency  of  data 
and  inconsistence  between  experiments.  General  conclusions  that  may  be  made 
on  the  microscale  effects  in  single  phase  flow  and  heat  transfer  include  [1]:  (1) 
friction  factors  in  laminar  microchannel  flow  are  apparently  lower  than  macroscale 
predictions  at  the  same  Re,  (2)  turbulent  heat  transfer  coefficients  are  apparently 
higher  than  theory  predicts  at  the  same  Re,  (3)  microscale  effects  seem  to  become 
more  prominent  for  channels  smaller  than  40/um,  (4)  slip  flow  theories  appears  to 
be  applicable  to  only  the  smallest  channels  in  gas  flow,  and  (5)  based  on  (1)  and 
(2),  the  Reynolds  analogy  appears  to  be  invalid  in  the  microscale  regime. 

Two-phase  flow  and  heat  transfer  with/without  phase  change  is  another  type  on 
microchannel  process  important  for  the  development  of  complete  systems.  The  ma¬ 
jor  advantage  of  two-phase  heat  transfer  with  phase  change  is  that  the  temperature 
gradient  along  the  length  of  flow  is  much  smaller  than  in  single-phase  flow.  This 
results  in  larger  values  of  the  heat  transfer  coefficient  and  reduces  the  requirement  of 
large  flow  rates.  Research  of  boiling  heat  transfer  in  open  and  closed  microchannels 
has  achieved  some  significant  results.  The  minimum  channel  size  for  which  data  are 
available  is  approximately  100/itm.  Maximum  heat  fluxes  of  350VF/cm^,  ISOW/cm^, 
and  SOW/crn^  have  been  found  for  open  and  closed  microchannel  convection  flow 
boiling,  and  micro  heat  pipes,  respectively  [3].  The  magnitude  of  these  heat  flux 
data  indicate  why  two-phase  heat  transfer  is  such  an  important  area  of  research. 
Some  of  the  most  promising  areas  for  future  development  are  the  areas  of  liquid 
metal  heat  pipes  and  micro  vapor  compression  cycles.  Convective  flow  boiling  in 
microgeometries  appears  to  be  a  technology  that  can  satisfy  the  demand  for  the 
dissipation  of  high  heat  fluxes  associated  with  electronic  and  laser  devices.  However, 
at  this  time,  the  amount  of  work  available  in  the  literature  related  to  open  and 
closed  microchannels  is  quite  limited.  Some  progress  concerning  convective  flow 
boiling  in  microgeometries  has  been  made  by  the  work  of  Zhukov  and  Yarmak  [21], 
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that  of  Peng  and  Wang  [2],  and  that  of  Bowers  and  Mudawar  [22].  Although  some 
obvious  differences  of  performance  compared  to  that  in  macrogeometries  have  been 
observed  in  their  work,  no  general  conclusions  can  be  drawn  in  this  area. 

Presently,  the  area  of  microscale  two  phase  flow  and  heat  transfer  without  phase 
change  is  being  dealt  with.  The  only  work  found  in  this  area  is  that  of  Suo  and 
Griffith  [23],  and  that  of  Stanley  [24]  . 

1.2  Objectives  of  Current  Resecurch 

As  mentioned  earlier,  work  in  microscale  condensation  flow  and  heat  transfer 
has  not  been  seen  in  the  literature.  As  a  result,  the  current  research  should  be 
expected  to  shed  considerable  light  on  this  field.  The  goal  of  this  research  is  to 
experimentally  investigate  condensation  flow  and  heat  transfer  of  steam  forced 
convection  in  a  horizontal  microtube.  A  stainless  steel  microtube  with  inside 
diameter  of  250/im  will  be  used  for  the  research.  The  experiments  will  be  expected 
to  determine  condensation  Nusselt  numbers  and  condensation  pressure  drops  in 
the  range  of  Reynolds  numbers  of  inlet  saturated  vapor  approximately  from  500 
to  5,000.  This  selected  range  of  Reynolds  numbers  is  due  to  the  consideration 
that  the  resulting  pressure  drops  would  be  appropriate  for  actual  applications. 
Experimental  outcomes  will  be  compared  to  macroscale  results  available  in  the 
literature  to  identify  possible  microscale  effects. 


CHAPTER  2 


REVIEW  OF  MACROSCALE 
CONDENSATION 

As  mentioned  in  Chapter  1,  there  have  not  been  any  experiments  reported  in  the 
literature  for  condensation  flow  and  heat  transfer  in  microgeometries.  This  chapter 
will  review  some  widely  used  and  accepted  results  acquired  for  macroscale  conden¬ 
sation,  which  will  be  used  for  comparison  with  the  current  experiment.  Types  of 
condensation  will  first  be  introduced,  and  then  correlations  for  condensation  heat 
transfer  and  the  method  to  obtain  condensation  pressure  drop  will  be  described. 

2.1  Types  of  Condensation 

Condensation  is  one  of  the  most  common  modes  of  heat  transfer.  Condensation 
is  defined  as  the  removal  of  heat  from  a  system  in  such  a  manner  that  vapor  is 
converted  into  liquid.  There  are  various  ways  in  which  condensate  may  be  formed. 
The  different  mechanisms  are  known  as  the  modes  of  condensation.  The  major 
modes  include  dropwise  condensation  and  film  condensation.  Since  it  is  impossible 
to  maintain  cooling  surfaces  in  a  state  that  maintains  dropwise  condensation  for 
long  periods  of  time,  film  condensation  is  of  practical  importance  in  applications 
for  industrial  condensers.  In  this  mode,  the  condensate  wets  the  surface  and  forms 
a  continuous  liquid  film. 

Film  condensation  can  occur  on  an  outside  surface  or  on  an  inside  surface  of 
a  tube.  For  the  latter  form,  film  condensation  takes  place  when  the  flowing  vapor 
is  driven  through  the  subcooled  inside  surface  of  the  tube.  Depending  on  the 
orientation  of  the  tube,  this  flowing  condensation  can  be  divided  into  horizontal, 
vertical,  and  inclined  condensation  in  the  tube.  Film  condensation  inside  horizontal 
tubes  will  be  of  interest  in  the  current  research. 
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2.2  Pattern  Identification  of 
Film  Condensation 

Flow  patterns  of  film  condensation  inside  a  horizontal  tube  depend  on  the 
distribution  of  forces  in  the  fluid  and  the  amount  of  condensate  in  a  cross  section. 
These  flow  patterns  are  shown  in  Figure  2.1  [25]. 

If  the  vapor  shear  stress  on  the  condensate  film  is  sufficiently  large  in  comparison 
to  the  effect  of  gravity,  and  if  the  amount  of  condensate  is  small,  then  there  will  be 
an  annular  flow. 

If  the  velocity  of  the  vapor  flow  is  moderate,  the  influence  of  gravity  will  cause 
the  condensate  to  flow  off  around  the  circumference  of  the  wall,  which  results  in  a 
condensation  pattern  called  stratified  flow. 

Spray  flow  is  the  result  of  shear  of  high  velocity  vapor.  In  this  pattern,  a  small 
amount  of  liquid  from  the  inside  surface  of  the  tube  is  entrained  into  the  core  of 
the  tube  and  forms  a  smog  mixture  of  vapor  and  liquid. 

Other  flow  patterns  include  intermittent  flow  and  bubble  flow.  An  intermittent 
flow,  which  consist  of  liquid  with  trapped  vapor  bubbles  or  of  a  liquid  slug  with 
vapor  zones  extending  to  the  wall,  forms  a  low  vapor  velocity  but  with  a  large  liquid 
component.  At  a  high  vapor  velocity  and  in  a  large  liquid  component,  a  bubble 
flow  with  bubbles  being  distributed  in  a  cohesive  liquid,  will  result. 

For  the  calculation  of  heat  transfer,  one  must  first  estimate  which  flow  pattern 
is  present.  A  starting  point  for  this  is  to  use  the  flow  pattern  chart  shown  in 
Figure  2.1.  In  this  chart,  the  term  (1  —  Q;)/Qf  in  the  horizontal  axis  represents  liquid 
volumetric  fraction  where  a  is  void  fraction  of  the  vapor  which  is  related  to  vapor 
quality  as 


Jg  the, vertical  axis  is  the  dimensionless  vapor  velocity,  which  is  the  ratio  of 

vapor  shear  axial  gradient  to  gravity  radial  gradient,  that  is 

•*  _ 

Aolgdipgipj  - 


(2.2) 


Dimensionless  Vapor  Velocity  j 


Figure  2.1.  Flow  Patterns  of  Condensation  Inside  Horizontal  Tubes 


9 


The  following  reference  values  can  be  used  to  differentiate  the  types  of  flow  [25]; 

Stratified  Flow:  <  1,  (1  —  a)/a  <  0.5 

Annular  Flow:  1  <  Jq  <  6,(1  —  a)/a  <  0.5 
Spray  Flow:  >  6,  (1  —  a)/Q;  <  0.5 

Slug  Flow:  Jq  <  0.01,  (1  —  «)/«  >  0.5 
Semi-Annular  Flow:  0.01  <  Jq  <  0.5,  (1  —  a)/a  >  0.5 
Bubble  Flow:  >  0.5,  (1  —  a)/a  >  0.5 

The  most  frequently  occurring  types  of  flow  are  stratified  flow  and  annular  flow. 
A  large  number  of  studies  available  in  the  literature  concern  these  two  kinds  of 
condensation  and  their  results  will  be  presented  in  the  following  sections.  It  should 
be  mentioned  that  spray  flow  condensation  would  definitely  have  a  higher  heat 
transfer  coefficient  than  both  stratified  flow  and  annular  flow  condensation  though 
no  data  for  spray  flow  have  been  found  in  the  literature. 


2.3  Stratified  Flow  Condensation 

Nusselt  first  proposed  (in  1916)  a  simple  theory  for  analyzing  the  laminar  film 
condensation  on  a  vertical  flat  surface  [26].  The  major  assumptions  made  by  Nusselt 
are  that  the  film  flow  is  laminar,  the  vapor  is  stagnant  and  exerts  no  drag  on  the 
downward  motion  of  the  condensate,  and  the  heat  transfer  through  the  film  is  by 
conduction  only.  By  solving  the  momentum  equation  separately  and  relating  the 
solutions  to  a  heat  balance,  Nusselt  obtained  the  average  Nusselt  number  as 


Nu  =  0.943 


'9P/(pf  -  Pa)kfh/g‘ 

U/J 

pjLAT 

T1/4 


(2.3) 


Nusselt ’s  film  condensation  theory  has  been  found  to  under-predict  the  actual 
Nu  in  some  situations.  The  discrepancy  has  been  primarily  attributed  to  the 
presence  of  waves  on  the  film  surface  at  low  film  Reynolds  numbers,  and  the  vapor 
shear  stress  at  the  vapor-condensate  film  interface  at  high  film  Reynolds  numbers 
[27]. 
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Condensation  occurring  inside  a  tube  differs  from  that  on  a  vertical  flat  surface 
in  that  the  vapor  entering  the  tube  with  some  magnitude  of  velocity  can  exert  a 
shear  stress  on  the  condensate  film.  For  low  vapor  velocity,  Chato’s  recommen¬ 
dation  [28]  for  the  average  condensation  Nusselt  number  inside  a  horizontal  tube 
is 


where 


Nu  =  0.555 


("-"l 

'gPfiPj  -  P9)kfh'jg 

U/j 

fifDAT 

1/4 


=  h,,  +  (3/8)c,,AT 


(2.4) 


is  the  revised  enthalpy  of  vaporization.  Chato’s  equation  is  limited  to  the  Reynolds 
numbers  of  the  inlet  saturated  vapor  in  the  range  Re^^i  <  35,000,  where  Rcg^i  is 
defined  as 


Re 

f^3,i 


(2.5) 


Similarity  can  been  seen  between  Chato’s  equation,  Eq.  (2.4),  and  Nusselt’s 
equation,  Eq.  (2.3).  An  ideal  stratified  flow  looks  like  that  shown  in  Figure  2.2. 
Due  to  the  low  vapor  velocity,  gravity  makes  the  condensate  formed  on  the  top 
of  the  inside  surface  of  the  tube  flow  off  around  the  circumference  downward  the 
bottom.  This  process  resembles  that  taking  place  on  a  vertical  flat  surface  analyzed 
by  Nusselt,  except  for  the  curvature  of  the  inside  surface  of  the  tube.  The  different 
coefficients  appearing  in  both  equations,  Eq.  (2.3)  and  Eq.  (2.4),  represent  this 
difference. 

It  should  be  noticed  that  stratified  flow  condensation  does  not  depend  on  Reynolds 
number  of  the  inlet  saturated  vapor  or  vapor  velocity  as  seen  in  Eq.  (2.4). 

By  far,  Chato’s  equation,  Eq.  (2.4),  is  the  most  widely  used  correlation  for 
stratified  flow  condensation. 


2.4  Annular  Flow  Condensation 

Another  common  type  of  film  flow  condensation,  annular  condensation,  has 
been  studied  by  a  number  of  researchers.  The  transition  from  stratified  flow  to 
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Figure  2.2.  Stratified  Film  Condensation 

annular  flow  is  caused  by  the  large  vapor  shear  stress.  There  are  five  widely-used 
empirical  or  semiempirical  correlations  available  in  the  literature  for  annular  flow 
condensation. 

Akers  et  al.  [29,  30]  used  four  dimensionless  parameters  to  correlate  the  average 
Nusselt  number  as: 


Nu  =  f  (Fry,  Reg,^q,  AT*,  Ref) , 


where  Prj  is  the  condensate  Prandtl  number, 


u,D  (p,\''^ 


Re  =  zm  ZL 
\Pa 

is  the  equivalent  vapor  Reynolds  number. 


AT*  =  —il— 

CgfAT 

is  the  dimensionless  subcooled  temperature,  and 


(2.6) 


(2.7) 


(2.8) 


Rej  =  ^  (2.9) 

is  the  film  Reynolds  number.  If  the  film  flow  is  laminar  (Re/  <  5,000),  Re/  is 
absent  in  Eq.  (2.6);  otherwise  all  four  parameters  should  be  employed. 

Akers  et  al.  [29,  30]  obtained  the  following  experimental  correlation: 
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For  laminar  flow  of  the  condensation  film,  that  is  Re/  <  5,000,  when  1,000  < 
R^g,eq  ^  20, 000,  the  average  Nusselt  number  is 

m  =  (2.10) 

and  when  20,000  <  Rcg^eq  <  100,000, 

WH  =  0.1Pry®AT*^/®i2e2g.  (2.11) 

For  turbulent  flow  of  the  condensation  film,  that  is  Re/  >  5, 000, 


Nu  =  0.026Pry"(Pe,,e,  +  Pe/)®-«,  (2.12) 

when  Rcg^eg  >  20, 000. 

The  second  correlation  for  the  local  Nusselt  number  of  annular  film  condensation 
is  that  of  Soliman  et  al.  [31]: 


p_0.65  n.i/2 

Nu  =  0.036— 

N 

where  Fo  is  the  wall  shear  stress  that  is  defined  as  F),  =  F/  +  where: 


(2.13) 


p  /  \  0.523  /  s  0.261' 

{  ri4  =  0-045Fe;«-2  x»-8°  +  5.70  ^  (1  -  ^ 

SmyTT^PgD*  \PgJ  ^  ’  \pjj 


+0.045Fe-°'2 

4f>*’ 


U1  ( 


0.105 


(1  -  ^ 

\pf. 


0.522' 


and 


Sm? /ir^pgD^ 


/D\  \  /  n  \  ^  /  \ 

=  -0.50  _ )  (2x  -  1  -  ^x)  ^  +  (i  _  3  +  2x)  ^ 

\LJ  \pfj  X  \pfj 


PI 


-0.50  (I)  [2(1  -  X)  +  (2^  -  ^  _  /,x)  (^)  +  2(1  -  X  -  /?  +  /3x) 

where  /?  is  a  constant  relating  interfacial  liquid  velocity  to  the  average  velocity 
throughout  the  condensate  film,  and 
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Re,,  =  ^  (2.14) 

y-g 

is  the  Reynolds  number  based  on  total  flow  rate  and  local  vapor  density. 

The  third  correlation  for  the  average  Nusselt  number  of  annular  film  condensa¬ 
tion  of  steam  was  proposed  by  Boyko  and  Kruzhilin  [32]: 


where 


Nu  =  0.024Re5;f  Pr?'*® 


2 


and 


(2.15) 


Res,t  = 


UffiD 


(2.16) 


is  the  Reynolds  number  assuming  all  mass  would  flow  as  liquid. 

Traviss  et  al.  [33]  reported  the  following  correlation  for  the  local  Nusselt  number: 


Nu  = 


/^o.g 
0.1  ^ 


-0.9 1 


1  -h  2.85 


where 


P2  =  0.707Pr/Pe?;^, 


for  Pe/,e,  <  50, 


P2  =  5Pr/  -I-  5/n[l  Pr/(0.09636Pe;;®f  -  1)] 


for  50  <  Pe/,e?  <  1125, 
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F2  =  5Prf  +  5ln{l  +  5Pr/)  +  2.5/n(0.00313Pe5;«f ) 
for  Rej,eq  >  1125,  and  Ref,eq  is  the  equivalent  Reynolds  number  defined  as 


Rej,eq  =  ^^(1  -  x).  (2.18) 

In  addition  to  the  above  four  correlations,  Shah  [34]  proposed  a  simple  equation 
which  is  considered  to  cover  the  widest  range  of  experimental  data.  This  equation 
represents  measurements  with  water,  i?  —  11,  R  —  12,  R  —  113,  methanol,  ethanol, 
benzene,  toluene,  and  trichloroethylene  during  condensation  in  vertical,  horizontal, 
and  inclined  tubes  with  inside  diameters  of  7.4  to  40mm.  The  range  of  experimental 
parameters  includes  the  reduced  pressure  0.002  <  Pr  <  0.44,  saturated  temperature 
21°C'  <  Tsat  <  OlO^C,  vapor  velocity  Zm/s  <  Ug  <  300m/s,  mass  flowrate 

lO.Skg/s  —  m^  <G  <  210.6A:5r/s  —  m^,  heat  flux  158iy/m^  <  9  <  1.893  x  10® 

Reynolds  number  100  <  Rej^t  <  63,000,  Prandtl  number  1  <  Pr  <  13,  and  vapor 
quality  0%  <  x  <  100%.  The  mean  error  of  this  equation  is  ±15.4%.  The  following 
is  the  Shah’s  equation  for  the  average  Nusselt  number: 


_  0.023Pe®/fPrS'' 

Nu  = - — i- 


X2  —  Xi 


(l  -  a:)'-®_^  3.8  ( 


1.8 


pO.38 


1.76 


0.04x2™ 

2.76  )  ’ 

/  J  a?! 


(2.19) 


where  pr  is  the  reduced  pressure  of  vapor,  and  Rej^t  is  the  Reynolds  number  defined 
in  Eq.  (2.16). 


2.5  Two  Phase  Pressure  Drop 

Two-phase  pressure  drop  in  a  horizontal  tube  is  the  result  of  the  action  of  two 
factors:  the  friction  force  and  momentum  changes.  The  general  equation  is 

AP  =  AP/  +  AP,„.  (2.20) 

The  pressure  drop  due  to  momentum  changes  is  given  by 


15 


w? 

^Pm  =  -  Vi).  (2.21) 

At  pressures  close  to  one  atmosphere  for  a  very  wide  variety  of  liquids,  the 
Lockhart  and  Martinelli  method  [35,  36]  is  considered  to  be  the  best  in  calculating 
two  phase  friction  pressure  drop.  The  friction  pressure  drop  is  determined  by  first 
calculating  the  pressure  drop  for  either  the  liquid  flowing  alone,  or  the  gas  flowing 
alone  in  the  tube  at  the  single  species  mass  flowrate.  That  is 


=  f(j 

Tng\ 

5)  2Al  ’ 

(2.22) 

or 

AP„  =  /(^ 

A  nif^Vf 

V  2Ag  ■ 

(2.23) 

In  the  above  two  equations,  /  is  the  friction  factor  that 

is  determined  from  the 

Moody  chart.  The  two-phase  friction  pressure  drop  is  then  obtained  by  multi¬ 
plying  either  the  liquid-only  pressure  drop,  or  the  gas-only  pressure  drop,  by  the 
appropriate  multiplier  or  <^f.  That  is 

AP/  =  <^jAPy,„  (2.24) 

or 

APf  =  (2.25) 

The  multiplier  (^^  or  <^y)  is  a  function  of  X  that  depends  on  the  flow  mechanism  of 
both  the  liquid  flowing  alone  and  the  gas  flowing  alone.  <f>g  or  (j)/  can  be  obtained 
either  from  the  Lockhart  and  Martinelli  chart  or  from  the  Lockhart  and  Martinelli 
data  [35,  36]  once  X  is  determined  from 

Re  J  Cg  y  Trig  j  Vg 

where  coefficients  m,  n,  C/,  and  Cg  are  determined  from  Table  2.1. 


(2.26) 
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Table  2.1.  Lockhart  and  Martinelli  Constants 


Constant 

Laminar  Regime 
Reg{orRef)  <  2,500 

Turbulent  Regime 

2,500  <  Reg{orRef)  <  50,000 

64 

0.316 

C'. 

64 

0.316 

m 

1 

0.25 

n 

1 

0.25 

2.6  Summary 

The  above  sections  have  offered  a  total  of  six  Nu  correlations  available  in  the 
literature  for  calculating  film  condensation  heat  transfer  inside  macrotubes,  one  for 
stratified  flow,  and  five  for  annular  flow.  They  are  seen  to  be  very  cumbersome  to 
use,  and  moreover  they  are  generally  coupled  with  pressure  drops.  Therefore,  some 
useful  points  are  summarized  here  for  the  research  of  macrotube  condensation  as 
follows: 

1.  The  minimum  inside  diameter  of  macrotubes  used  to  acquire  those  six  heat 
transfer  correlations  is  7.4mm. 

2.  There  are  up  to  six  types  of  Reynolds  numbers  defined  and  used  (Eq.  (2.5), 
Eq.  (2.7),  Eq.  (2.9),  Eq.  (2.14),  Eq.  (2.16),  and  Eq.  (2.18))  to  relate  Nusselt 
numbers.  Since  these  different  types  of  Reynolds  number  are  usually  not 
compatible  with  each  other,  the  different  results  of  Nusselt  numbers  reported 
are  difficult  to  compare. 

3.  Eight  dimensionless  parameters  are  used  to  relate  Nusselt  numbers.  These 
Nusselt  numbers  may  be  expressed  symbolically  as 

Nuor  Nu  =  /{six  types  of  Re,x,AT*,Pr,d/L,Pr,pglpf,pgfpf).  (2.27) 

4.  The  correlation  for  stratified  flow  condensation  (Eq.  (2.4))  is  independent 
of  Reynolds  numbers  or  velocities  of  the  inlet  saturated  vapor.  The  five 
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correlations  for  annular  flow  condensation  (Eq.  (2.10-2.12),  Eq.  (2.13),  Eq. 
(2.15),  Eq.  (2.17),  and  Eq.  (2.19))  all  depend  on  Reynolds  number. 

5.  Spray  flow  condensation  should  have  a  higher  heat  transfer  coefficient  than 
both  stratifled  and  annular  condensation  though  no  correlations  for  this  type 
of  condensation  are  available  in  the  literature. 

6.  The  vapor  quality  is  present  in  all  the  five  correlations  for  annular  flow 
condensation  (Eq.  (2.10-2.12),  Eq.  (2.13),  Eq.  (2.15),  Eq.  (2.17),  and  Eq. 
(2.19)). 

7.  The  heat  transfer  and  pressure  drops  are  generally  coupled. 

8.  Due  to  the  change  of  vapor  quality  along  the  length  of  a  tube,  a  numerical 
integration  must  generally  be  executed  between  the  inlet  and  outlet  of  the 
tube  to  acquire  the  overall  two-phase  pressure  drop  and  the  average  Nusselt 
number  from  the  correlations  for  the  local  Nusselt  number.  That  is,  the  tube 
is  discretized  into  a  series  of  short  sections  and  the  separate  contributions  of 
pressure  drop  and  heat  transfer  of  each  section  are  summed. 

9.  Errors  for  all  those  heat  transfer  correlations  are  reported  to  be  quite  large. 
A  100%  or  even  larger  error  is  a  common  situation  [36,  37]. 

10.  The  effect  of  Reynolds  number  on  Nusselt  number  is  observable  only  for  a  large 
Re  range.  For  instance,  the  Re  range  in  Akers  et  al.  [29,  30]  is  10^  <  Re  < 
10®.  This  result  is  due  to  the  large  uncertainty  associated  with  condensation 
experiments  which  tends  to  mask  the  Re  effect  when  only  a  small  Re  range 
is  considered.  For  microscale  tube  condensation,  the  Re  range  is  expected  to 
be  small  (i.e.,  500  <  Re  <  5,000)  due  to  pressure  drop  constraints  that  may 
make  it  difficult  to  observe  the  Re  effect  on  Nu. 


CHAPTER  3 


EXPERIMENTS  AND  MODELING 

The  objective  of  the  current  research  is  to  experimentally  investigate  the  heat 
transfer  and  pressure  drop  of  steam  condensation  inside  a  horizontal  microtube. 
This  chapter  will  describe  the  experimental  apparatus  and  procedure.  Then,  cor¬ 
relations  of  the  data  for  average  Nusselt  number  and  pressure  drop  will  be  derived. 

3.1  Experimental  Apparatus 

The  experimental  apparatus  is  shown  schematically  in  Figure  3.1,  and  a  photo¬ 
graph  of  this  experimental  apparatus  is  shown  in  Figure  3.2.  There  are  two  loops 
in  the  apparatus.  The  first  flow  loop  is  used  to  supply  water  (from  the  building 
supply)  that  is  to  be  condensed  after  it  is  heated  to  superheated  steam.  The  second 
loop  is  used  to  supply  the  cooling  water,  also  from  tap  water,  that  is  used  to  cool 
the  microtube  in  which  condensation  is  taking  place.  All  tubes  in  the  two  loops 
are  stainless  steel  with  6.35mm  (0.25in)  outside  diameter.  The  two  streams  of  the 
tap  water  are  filtered  by  two  7fj,m  filters,  respectively,  before  the  water  enters  each 
loop.  The  microtube  has  inside  diameter  of  250/xm,  outside  diameter  of 
and  length  of  4877im.  Table  3.1  lists  the  main  performance  data  of  instruments 
used  in  the  experiment. 

The  superheat  temperature  of  steam  at  the  inlet  of  the  microtube  is  controlled 
by  a  combination  of  a  temperature  controller  and  pressure  transducer.  The  tem¬ 
perature  controller  is  connected  to  a  thermocouple  that  measures  the  temperature 
of  steam  at  the  inlet  of  the  microtube.  This  temperature  is  fed  back  to  control 
the  on-off  of  a  rope  heater.  The  pressure  transducer  measures  gauge  pressure 
and  thus  the  absolute  pressure  at  the  the  inlet  of  the  microtube  is  obtained  by 
adding  the  inlet  gauge  pressure  to  the  outlet  absolute  pressure  that  is  equal  to  the 
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Pressure  Gauge 


Pressure  Gauge 


Figure  3.1.  Schematic  of  Experimental  Apparatus 
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Figure  3.2.  Photograph  of  Experimental  Apparatus 
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Table  3.1.  Limiting  Conditions  of  Main  Equipment 


Component  Name, 
Type  and  Maker 

Pressure  PS  I 
(kPa) 

Temperature 

(°C) 

Flow  Rate 
(ml/min) 

Accuracy 

Pressure  Transducer 
PX202-015GV 
OMEGA 

15(103) 

-40  ~  125‘’C' 

N/A 

±0.25% 

Pressure  Indicator 
DP41-S 

OMEGA 

N/A 

N/A 

N/A 

l^V 

Thermocouple 

5TC-TT-J-36-36 

OMEGA 

N/A 

(Type  J) 

N/A 

±0.1%reading 

+0.6°C 

Thermal  Meter 

HH23 

OMEGA 

N/A 

N/A 

N/A 

dtO. iVoreading 
-f0.6®C  (including 
thermocouple) 

Temperature  Controller 
CN3240 

OMEGA 

N/A 

-73  760®  a 

N/A 

±0.14%o /span 
+  1  digital 

Flow  Meter 
(Condensing  Loop) 
GF7060  GILMONT 

125  (862) 

N/A 

0.0091  ~  1.28(G/ass) 
0.0576  ~  5.49(S.S.) 

±2%reading 

Flow  Meter 
(Cooling  Loop) 
GF7160  GILMONT 

125  (862) 

N/A 

0.0212  ~  4.50(Giass) 
0.115  ~  16.6(S.S.) 

±2%readmp 

Pressure  gauge 
PGC-20L-100 
OMEGA 

100(690) 

N/A 

N/A 

N/A 

atmospheric  pressure  in  the  test  field  (the  outlet  of  the  microtube  is  open  to  the 
atmosphere).  The  atmosphere  pressure  is  determined  by  using  a  barometer  placed 
in  the  test  field.  Volumetric  flowrates  of  tap  water  for  both  loops  are  measured 
by  two  fine  float  flowmeters,  respectively.  Figure  3.3,  Figure  3.4  and  Figure  3.5 
depict  the  heat  exchanger  in  which  condensation  occurs.  The  heat  exchanger  is  a 
non-mixing  cross-flow  heat  exchanger.  The  outside  shell  is  a  Swagelok  cross  fitting, 
and  the  inside  tube  is  the  microtube  being  studied.  The  cooling  water  flows  over 
the  microtube.  The  amount  of  condensed  heat  removed  is  calculated  using  the 
change  of  enthalpy  of  the  cooling  water.  Two  thermocouples  located  at  the  inlet 
and  outlet  of  the  heat  exchanger,  respectively,  measure  temperatures  that  are  used 
to  calculate  the  change  of  enthalpy  of  the  cooling  water.  The  outside  surface  of 
the  heat  exchanger  is  well-insulated  by  the  wrapped  insulation  tape,  thus  heat  loss 
through  the  heat  exchanger  may  be  ignored. 

Outside  wall  temperatures  of  the  microtube  are  measured  by  four  fine  ther- 
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Cooling  Warer  In 


Figure  3.3.  Schematic  of  Heat  Exchanger 

mocouples  that  are  glued  to  the  microtube  outside  wall  and  are  equally  spaced 
along  the  flow  direction.  The  significant  error  of  measurement  of  the  outside  wall 
temperatures  may  be  ignored  because  the  bead  diameter  of  those  fine  thermocou¬ 
ples  {127 fim)  is  small  enough  compared  to  the  outside  diameter  of  the  microtube 
(790/xm).  All  thermocouples  used  in  the  experiment  are  connected  to  a  thermal 
meter. 


3.2  Experimental  Procedure 

The  experiments  are  conducted  according  to  the  following  procedure; 

1.  Evacuate  air  in  the  condensing  loop  by  letting  the  tap  water  flow  through  the 
loop  at  a  large  flow  rate  for  a  long  period  of  time  (approximately  5  minutes). 
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Figure  3.4.  Photograph  of  a  Portion  of  Experimental  Apparatus 


Figure  3.5.  Photograph  of  Heat  Exchanger 


24 


2.  Set  the  volumetric  flowrate  into  the  condensing  loop  at  a  fixed  value  such  that 
an  expected  Reynolds  number  (Re^ ,)  of  inlet  saturated  vapor  is  achieved.  The 
expected  range  of  Rcg  i  is  approximately  from  500  to  5, 000. 

3.  Regulate  the  volumetric  flowrate  into  the  cooling  loop  such  that  an  expected 
vapor  quality  (2:2)  at  the  outlet  of  the  microtube  is  achieved  at  the  fixed  Rtg^i. 
The  expected  range  of  X2  is  approximately  from  0.1  to  0.9. 

4.  Obtain  an  expected  superheat  temperature  of  the  inlet  vapor  by  setting 
a  value  of  the  temperature  controller,  in  comparison  with  the  reading  of 
the  pressure  transducer.  The  superheat  temperature  of  the  inlet  vapor  is 
controlled  to  be  less  than  2b°C. 

5.  Collect  data  when  a  steady  state  has  been  achieved.  Collected  data  include 
those  which  will  be  shown  in  next  section. 

6.  Change  the  volumetric  flowrate  into  the  condensing  loop  to  another  value 
which  corresponds  to  another  Rcg  i. 

7.  Repeat  steps  2  ~  6  until  all  the  values  of  the  expected  Rtg^i  and  X2  are 
considered. 

3.3  Experimental  Modeling 

3.3.1  Directly  Measurable  Quantities  in  the  Experiment 

Referring  to  Figure  3.1  and  Figure  3.3,  15  parameters  are  directly  measurable 
by  means  of  the  data  acquisition  instruments.  These  15  parameters  are  shown  in 
Table  3.2. 


3.3.2  Inlet  Vapor  Quality  and  Thermodynamic 

Process  of  Steam 

The  steam  is  heated  to  a  superheated  state  to  ensure  that  only  vapor  is  intro¬ 
duced  to  the  microtube  where  condensation  occurs.  Referring  to  Figure  3.3,  an 
assumption  is  made  that  the  wall  between  the  supply  tubing  and  the  microtube 
is  ideally  thermally  insulated,  then  the  superheated  steam  will,  when  condensed 
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Table  3.2.  Directly  Measurable  Quantities  in  the  Experiment 


DEB 

Variable 

Description 

1 

inside  diameter  of  the  microtube 

mm 

do 

outside  diameter  of  the  microtube 

IBB 

L 

length  of  the  microtube 

IDl 

inside  diameter  of  the  supply  tubing 

5 

IBHH 

inlet  temperature  of  water  into  both  loops 

6 

cooling  water  temperature  at  the  outlet  of  the  heat  exchanger 

Tspht 

superheat  temperature  of  steam 

Qi 

mmmm 

wall  temperature  of  the  microtube  at  point  1  (Figure  3.2) 

BB 

TwI,o,2 

wall  temperature  of  the  microtube  at  point  2  (Figure  3.2) 

10 

wall  temperature  of  the  microtube  at  point  3  (Figure  3.2) 

11 

EiUM 

wall  temperature  of  the  microtube  at  point  4  (Figure  3.2) 

wm 

Po 

atmospheric  pressure 

P spht 

gauge  pressure  of  steam  in  the  inlet  of  the  microtube 

14 

Vo, 

volumetric  flowrate  of  water  entering  the  condensation  loop 

15 

volumetric  flowrate  of  water  entering  the  cooling  loop 

inside  the  microtube,  arrive  at  a  saturated  vapor  state  somewhere  downstream 
close  to  the  inlet  of  the  microtube.  Let  us  call  this  length  the  superheated  length, 
Lsphti  which  is  marked  in  Figure  3.3.  Due  to  the  good  thermal  conductivity  of 
the  stainless  steel  wall  between  the  supply  tubing  and  the  microtube,  there  would 
be  some  unavoidable  amount  of  heat  entering  the  cooling  water  through  the  wall. 
Accordingly,  the  superheated  length,  Lspht,  would  be  reduced  in  the  actual  situation. 
If  the  superheat  temperature  of  the  inlet  vapor  is  not  too  high,  Lgpht  would  be  very 
close  to  zero  and  the  location  where  the  saturated  vapor  state  of  steam  is  achieved 
would  be  very  close  to  the  inlet  of  the  microtube.  In  the  experiment,  the  superheat 
temperature  will  be  controlled  to  be  less  than  25“C'.  Under  this  limiting  value,  the 
following  calculation  will  show  that  Lspht  is  so  small  that  it  can  be  assumed  to  be 
zero. 

It  is  assumed  that: 

1.  the  wall  between  the  supply  tubing  and  the  microtube  is  ideally  thermally 
insulated. 
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2.  the  condensation  heat  removal  along  the  length  of  the  microtube  is  uniform, 

3.  the  absolute  pressure  across  the  microtube  is  O.lMPa, 

4.  the  superheat  temperature  of  steam  at  the  inlet  of  the  microtube  (point  1') 
is  25‘'C, 

5.  the  vapor  quality  (X2)  at  the  outlet  of  the  microtube  (point  2)  is  from  0.2  to 
0.8,  and 

6.  the  change  of  kinetic  energy  from  point  1'  to  point  1  is  ignored. 

With  these  conditions,  we  have 

Lspht  Q  spht 

^  Qt 

_  ^cd{hi 

mcd{hi  -  /i2) 

—  ~  hsat,l 

hi  —  /i2 

_  hi  hgdf^i 

{hi  hsat^l)  “h  {hsatjl  h2) 

_  hi  hsai,l 

{hi  —  hsat.l)  +  —  {hsatj2  +  ^2hsatjg2)] 

50 

~  50  +  [2675  -  (417  +  0.8  x  2258)] 

=  10.0%,  (3.1) 

where  values  of  enthalpy  of  steam  are  taken  from  [38],  and  X2  is  assumed  to  be  0.8 
which  would  be  the  worst  situation. 

It  is  noticed  that  the  above  calculation  is  based  on  the  ideal  insulation  as¬ 
sumption  (see  assumption  1).  However,  the  actual  process  with  nonideal  insulation 
would  reduce  the  superheated  length  to  a  smaller  value.  Therefore,  it  seems  to  be  a 
reasonable  assumption  that  the  steam  achieves  a  saturated  state  at  the  inlet  (point 
1)  of  the  microtube,  that  is,  the  inlet  vapor  quality  xi  w  100%. 

On  the  other  hand,  for  the  sake  of  accuracy,  the  amount  of  heat  required  for 
condensation  between  the  inlet  (point  1)  and  outlet  (point  2)  of  the  microtube. 
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which  will  be  used  to  calculate  average  condensation  Nusselt  numbers,  should  be 
obtained  by  deducting  the  heat  gain  through  the  wall  between  the  supply  tubing 
and  the  microtube  from  the  total  amount  of  heat  removed  by  the  cooling  water. 

The  process  from  point  1'  to  point  1  is  a  single  phase  process  in  which  the 
superheated  vapor  state  at  point  1'  changes  to  saturated  vapor  state  at  point  1. 

The  process  from  point  2  to  point  0  is  a  process  of  two  phase  flow'  without 
phase  change,  in  which  the  steam  jets  into  the  atmosphere  from  the  outlet  of  the 
microtube.  It  is  assumed  that  the  distance  between  point  2  and  point  0  is  so  short 
that  this  process  can  be  considered  to  be  thermally  insulated  from  the  atmosphere. 
It  is  also  reasonable  to  assume  that  the  cross-sectional  area  of  stream  of  the  steam 
at  point  0  is  the  same  as  that  of  the  microtube  (^40),  and  the  absolute  pressure  at 
point  0  is  equal  to  the  atmospheric  pressure. 

It  is  noticed  that  there  will  be  a  sudden  contraction  for  the  fluid  from  point  1' 
to  point  1  and  a  sudden  expansion  from  point  2  to  point  0.  Therefore,  the  overall 
thermodynamic  process  that  the  steam  undergoes  from  point  T  to  point  0  should 
be: 

1.  from  point  1'  to  point  1:  single  phase  heat  transfer  without  phase  change,  a 
sudden  contraction  of  flow  involved. 

2.  from  point  1  to  point  2:  condensation  phase  change. 

3.  from  point  2  to  point  0:  two  phase  adiabatic  flow  without  phase  change,  a 
sudden  expansion  of  flow  involved. 

The  following  is  the  derivation  of  data  reduction  equations  for  IVu,  X2,  and  AP. 
It  will  be  seen  from  these  equations  that  they  are  coupled. 

3.3.3  Pressure  Drop 

Applying  an  energy  balance  between  point  1'  and  point  1  (see  Figure  3.3)  gives 
an  equation  for  the  saturated  vapor  pressure  {Psat,i)  at  the  inlet  of  the  microtube: 


28 


Pspht  +  =  Psat,l  +  (1  +  (3.2) 

where  Kc  is  the  contraction  coefficient  due  to  the  sudden  contraction  from  diameter 
D  to  di.  Empirical  correlations  are  available  for  such  as  one  presented  by  White 

[39]: 


A'  =  0.42  1 


Based  on  mass  continuity  of  steam, 


~  Aovi  ’ 

and  the  velocity  of  superheated  steam,  Ui,  is  given  by 

mcdVi  . 

u.  =  — .  (3.5) 

The  saturated  vapor  pressure  at  the  inside  inlet  of  the  microtube  is  obtained  by 
inserting  Eqs.  (3.3),  (3.4),  and  (3.5)  into  Eq.  (3.2): 


1  frhcdViYl 


In  order  to  find  the  pressure  {Psat,2)  at  the  outlet  of  the  microtube,  the  vapor 
quality  at  the  outlet  of  the  microtube  must  be  determined  first.  By  applying  an 
energy  balance  between  point  T  and  point  0  (see  Figure  3.2),  we  have 

*  11^ 

rhcdihi  +  -^)  =  Qt  +  mcd{h,at,o  +  y),  (3.7) 

where  Qt  is  the  total  rate  of  heat  removed  by  the  cooling  water,  that  is 


Mass  continuity  of  the  steam  yields 


—  — 

U\  AqVi  ’ 

where  Usat.o  is  the  specific  volume  at  point  0  that  is  given  by 
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‘^sat,0  —  '^sat,fO  “H  ^0'^sat,Jg0^  (3.10) 

and  the  velocity,  Ui  in  Eq.  (3.7),  has  been  determined  by  Eq.  (3.5).  The  enthalpy 
at  point  0  is 

hsat,0  —  ^sat,  fO  ^0^sat,/g0-  (^■^^) 

By  substituting  Eqs.  (3.8)  to  Eq.  (3.11)  into  Eq.  (3.7),  Xq  can  be  determined 
from  the  following  equation: 

01X0  +  02X0  +  03  =  0,  (3.12) 

where  oi,  02,  and  03  are  three  constants: 


1  f  ^cd'^satffgO 


02  —  hsatJgO  +  (viVsatJoVsatJgo), 

^  k.  \  \  U  U  ^  1  f  '^cd'^satJo\ 

<'‘""-'‘"■■■>+'‘••'■'“-'‘>-21—1  +21—4^)  ■ 

The  equation  for  the  pressure  (Psata)  at  the  outlet  of  the  microtube  is  obtained 
by  applying  the  energy  balance  between  point  2  and  point  0: 


Psat,2  + 


»/2  n.2 

^sat,2  _  o  1  ^0  ,  ^sat,2 

n -  —  ^0  +  Z - h  l\e- - : 

^^50/, 2  ^Vsai,0  ^Vsat,2 


(3.13) 


where  is  the  expansion  coefficient  due  to  the  sudden  expansion  from  diameter 
di  to  the  infinite  space  (the  atmosphere),  which  is  determined  to  be  [39] 


A'e  =  1. 

The  velocity,  uq  in  Eq.  (3.13),  is  given  by 


(3.14) 
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Uo  = 


'^cd'^satyO 

Ao 


(3.15) 


By  considering  Eq.  (3.10),  and  combining  Eqs.  (3.14)  and  (3.15)  into  Eq.  (3.13), 
we  obtain 


Psai,2  —  ^0  + 


' fmcd\  / 

j  {VsfxtjQ  4“  ^O'^satJgO)  •  (3.16) 


2(Vsa<,/0  "b  ^Q'^sat,}go) 

The  condensation  pressure  drop  is  then  acquired  from  Eq.  (3.6)  and  Eq.  (3.16): 


^Pcd  —  Paat,\  —  Paat,2-  (3-17) 

3.3.4  Outlet  Vapor  Quality 

The  result  for  vapor  quality  {X2)  at  the  inside  outlet  of  the  microtube  is  obtained 
directly  from  an  energy  balance  between  point  2  and  point  0,  that  is 


^2 


haat,J0  ~b  ^phsatifgO  ^sa<,/2 
h'aai,fg2 


(3.18) 


It  is  noticed  that  both  hsat,f2  ^'Hd  hsatjg2  a  function  of  Paat,2i  which  has  been 
determined  in  Eq.  (3.16). 


3.3.5  Nusselt  Number 

Finally,  the  condensation  average  Nusselt  number  is  determined  from  the  fol¬ 
lowing  definition: 


Nu  = 


hdi 


(3.19) 


where  h  is  the  average  heat  transfer  coefficient  determined  from 


jr _  Qi  Qspht 

TvdiLAT  ’ 


(3.20) 


where  Qspht  is  the  rate  of  heat  into  the  cooling  water  from  the  wall  between  the 
supply  tubing  and  the  microtube  as  explained  previously.  Qspht  can  be  determined 
from  an  energy  balance  between  point  1'  and  point  1: 
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spht 


rhcd 

,  ,  uf 

^sat,l  H  ^ 

2 

\  Ui  J  \ 

(3.21) 


In  Eq.  (3.20),  AT  is  the  average  subcooled  temperature  given  by 


(3.22) 

where  the  inside  wall  mean  temperature  of  the  microtube  Tu;i,i  can  be  obtained 
from  the  measured  outside  wall  temperatures  by  considering  the  thermal  conduction 
between  the  inside  wall  and  outside  wall  of  the  microtube: 


=  2  +  (<3,  -  (3.23) 

^  j  =  l  L'KKyjlL 

The  average  Nusselt  number  is  determined  by  combining  Eqs.  (3.20),  (3.21), 

(3.22)  and  (3.23)  into  Eq.  (3.19). 


CHAPTER  4 


RESULTS  AND  DISCUSSION 

The  macroscale  empirical  and  semiempirical  correlations  that  are  most  widely 
used  for  film  condensation  in  horizontal  tubes  have  been  presented  in  Chapter  2. 
This  chapter  will  first  compare  these  correlations  with  the  results  obtained  in  the 
current  experiment.  After  a  comparison  is  made,  a  significant  deviation  between 
them  will  be  seen.  On  the  other  hand,  it  will  be  seen  that  conventional  methods 
for  Nusselt  number  correlations  seem  not  to  work  for  the  current  experimental 
results.  It  is  found  that  the  correlated  Nusselt  numbers  in  the  current  research  are 
dependent  on  one  parameter  only  instead  of  many  parameters  as  is  expected.  Some 
attempts  to  explain  these  significant  differences  will  be  presented  in  the  latter  part 
of  this  chapter. 

4.1  Ranges  of  Experimental  Pcirameters 

Forty  six  sets  of  data  have  been  collected  and  are  listed  in  Appendix  B.  Ranges 
of  experimental  parameters  are  shown  in  Table  4.1.  The  superheat  temperature  of 
steam  {Tspht)  at  the  inlet  of  the  microtube  is  on  the  order  of  25° (7,  indicating  that 
it  is  reasonable  to  employ  the  assumption  that  the  inlet  vapor  quality  (xj)  is  100% 
as  discussed  in  Chapter  3.  The  velocity  of  the  inlet  saturated  vapor  {ug^i)  is  up  to 
276m/s  that  may  cause  heating  effects.  This  will  be  discussed  later  in  this  chapter. 

4.2  Flow  Patterns  Based  on  the 
Macroscale  Method 

One  method  to  determine  the  flow  pattern  in  the  current  experiment  is  from  the 
flow  pattern  map  shown  in  Figure  2.1.  From  the  values  of  the  dimensionless  vapor 
velocity  and  liquid  component  (1  -  a) /or  in  Table  4.1,  it  is  seen  that  current 
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Table  4.1.  Ranges  of  Experimental  Parameters 


393  ~  4683 

Xi 

100(%) 

33  ~  276(m/s) 

X2 

24.3  ~  94.2(%) 

AT 

6.3  ~21.3(“C') 

3h 

14.8  ~  150.9 

Tspht 

12.2  ~  24.5(°C') 

(1  —  a)la 

2.Z  X  10-“*  ~  41  X  lO-'* 

AP 

-3.4  ~  65.3(A:Pa) 

experimental  results  should  be  located  in  the  regime  of  spray  flow.  On  the  other 
hand,  Chato  recommended  [28]  that  for  stratified  flow,  Reynolds  numbers  of  inlet 
saturated  vapor  be  <  35, 000.  From  this,  the  current  experimental  results 
would  most  likely  be  in  the  stratified  flow  regime.  Because  the  experimental  or 
theoretical  correlations  of  Nusselt  numbers  for  spray  flow  are  non-existent  in  the 
literature,  and  because  the  annular  flow  regime  is  located  between  the  stratified  flow 
and  spray  flow  regimes  as  shown  in  Figure  2.1,  it  should  be  reasonable  to  choose 
available  correlations  from  both  stratified  flow  and  annular  flow  as  references  for 
comparing  current  experimental  results.  Note  that  no  attempt  was  made  to  verify 
visually  the  type  of  flow  that  occurred  in  the  tests. 

4.3  Comparisons  of  Experimental  Results  with 
Macroscale  Correlations 

The  experimental  results  for  Nusselt  numbers,  pressure  drops,  and  outlet  vapor 
quality  are  acquired  based  on  the  models  established  in  the  Chapter  3.  Properties 
of  water  are  obtained  by  a  curve-fitting  program  that  takes  data  from  the  property 
tables  in  the  literature  [38,  40).  The  formulas  used  for  the  curve-fittings  are  polyno¬ 
mials.  The  order  of  each  polynomial  is  determined  based  on  (1)  the  range  of  data 
to  be  used,  (2)  the  nature  of  the  data  in  the  tables,  and  (3)  a  trial-and-error  basis. 
Properties  of  stainless  steel  of  the  microtube  being  studied  are  from  MicroGroup, 
Inc.  [41]. 


34 


4.3.1  Nusselt  Number 

Four  empirical  and  semiempirical  correlations  for  film  condensation  Nusselt 
numbers  have  been  chosen  to  make  comparison  with  current  experimental  results. 
Refer  to  Chapter  2  where  the  correlations  that  are  most  widely  used  in  the  literature 
have  been  presented.  The  chosen  correlations  for  this  comparison  include  those  of 
Chato  [28],  Akers  et  al.  [29,  30],  Boyko  et  al.  [32],  and  Shah  [34].  It  is  noticed 
that  the  correlation  of  Chato  is  suitable  for  film  condensation  in  stratified  flow 
in  horizontal  tubes,  and  the  other  three  correlations  chosen  here  are  for  annular 
flow  in  horizontal  or  nonhorizontal  tubes.  As  mentioned  above,  it  is  unclear  which 
types  of  film  condensation  occurred  in  the  current  experiment,  therefore  the  above 
correlations  which  cover  a  wide  range  from  stratified  flow  to  annular  flow  are 
accordingly  selected. 
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Figure  4.1.  Comparison  of  Average  Experimental  Nusselt  Number  with  Correla¬ 
tions  in  the  Literature 
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The  Nu  data  from  the  experiments  are  compared  to  the  macroscale  correlations 
in  Figure  4.1.  In  this  comparison,  the  Reynolds  number  of  the  inlet  saturated  vapor 
is  selected  as  a  common  coordinate.  It  is  seen  that  the  current  experimental  data 
are  much  less  than  all  the  macroscale  correlations.  Chato’s  correlation  appears  to 
be  the  best  of  all  the  correlations  in  predicting  the  experimental  data.  Upon  closer 
examination  of  the  Chato  correlation,  revealed  in  Figure  4.2,  the  departure  is  still 
very  apparent.  The  value  of  Nusselt  number  averaged  over  the  overall  range  of  the 
experimental  Reynolds  numbers  is  approximately  3.06.  The  Nusselt  numbers  of  the 
chosen  correlations  for  the  current  comparison  from  that  of  Chato,  Shah,  Boyko  et 
al.,  and  Akers  et  al.  are  approximately  3,  10,  20  and  50  times  higher,  respectively, 
than  the  averaged  experimental  Nusselt  number  (3.06). 


Figure  4.2.  Comparison  of  Experimental  Nusselt  Numbers  with  Chato’s  Correla¬ 
tion 
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Figure  4.3.  Comparison  of  Pressure  Drops 

4.3.2  Pressure  Drop 

A  comparison  of  experimental  pressure  drop  with  that  predicted  by  the  Lockhart 
and  Martinelli  method  [35]  is  shown  in  Figure  4.3.  It  is  seen  that  the  experimental 
pressure  data  are  mostly  located  in  the  range  of  prediction,  which  account  for  about 
85%  of  the  total  data. 

4.4  Correlation  of  Experimental  Nusselt  Numbers 

Conventional  data  correlating  methods  for  film  condensation  in  horizontal  tubes 
have  been  attempted  to  correlate  the  current  experimental  Nusselt  numbers,  but 
they  appear  to  not  be  applicable. 

First,  no  obvious  direct  relationship  between  the  average  Nusselt  numbers  (Nu) 
and  Reynolds  numbers  is  found,  and  then,  the  method  used  by  Akers  et  al. 

[29,  30]  is  also  attempted. 
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As  described  in  Chapter  2,  Akers  et  al.  [29,  30]  used  four  dimensionless  parame¬ 
ters  to  correlate  their  data  of  the  average  Nusselt  number  for  annular  flow  as  shown 
in  Eq.  (2.6).  When  the  film  flow  is  laminar  {Rej  <  5, 000)  and  the  Prandtl  number 
is  constant  as  is  the  case  in  the  current  experiment,  four  dimensionless  parameters 
used  by  Akers  et  al.  can  be  reduced  to  two.  That  is,  the  average  Nusselt  number 
has  the  the  form  of 

where  A,  B,  C  are  constants.  Unfortunately,  this  method  does  not  produce  a  good 
correlation  in  the  current  experiment. 

One  possible  explanation  for  the  above  results  is  that  the  range  of  Reynolds 
number  (approximately  500  <  <  5,000)  in  the  current  experiment  is  so 

narrow  that  no  obvious  trend  can  be  observed,  given  the  large  data  uncertainty. 
This  can  be  the  similar  situation  to  that  occurring  in  macroscale  condensation 
experiments  as  discussed  in  Chapter  2. 

The  only  obvious  relationship  found  for  the  experimental  Nusselt  numbers  that 
reasonably  correlates  the  data  takes  the  form  of 

Nu  =  9.60Ax*•“^  (4.1) 

where  Ax  is  the  change  of  vapor  quality  from  the  inlet  to  outlet  of  the  microtube. 
It  is  noted  that  this  correlation  is  almost  linear.  Fig  4.4  depicts  this  correlation.  It 
is  seen  that  the  experimental  data  are  located  in  the  range  of  -50%  to  -f  75%  of 
the  above  correlation. 

It  was  rhentioned  previously  that  all  the  correlations  available  in  the  literature 
and  presented  in  Chapter  2,  except  that  of  Chato,  include  vapor  quality  as  one  of 
dimensionless  variables  for  correlating  Nusselt  numbers  (see  Eq.  (2.27)).  Therefore, 
the  correlation  currently  obtained  for  Nusselt  numbers,  Eq.  (4.1),  could  be  used  as 
a  correlation  in  the  design  of  condensers.  The  sole  dependence  of  Nusselt  number  on 
the  change  of  vapor  quality  found  in  the  current  experiment  simplifies  its  application 
in  a  design  process.  On  the  other  hand,  the  very  small  range  of  parameters  tested 
limits  its  general  applicability. 
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Figure  4.4.  The  Effect  of  Quality  Change  on  Nusselt  Number 

4.5  Correlation  of  Dimensionless  Temperature 

vs  Reynolds  Number 

Experimental  data  reveal  that  dimensionless  temperature  can  be  related  to 
Reynolds  number  of  the  inlet  saturated  vapor  as 

AT*  =  2.17  X  (4.2) 

Figure  4.5.  shows  this  correlation.  Experimental  data  are  located  in  the  range 
of  —35%  to  +70%  of  the  correlation. 

As  discussed  previously,  no  dependence  of  Nu  on  Rtg,i  is  observed  in  the 
experiment.  Thus,  the  average  condensing  heat  transfer  coefficient  (h)  may  be 
assumed  to  be  a  constant  in  the  range  of  Rtg^i  tested.  On  the  other  hand,  it  was 
observed  in  the  experiment  that  the  heat  flux  {q)  removed  by  the  cooling  water 
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Figure  4.5.  The  Effect  of  Inlet  Saturated  Vapor  Reynolds  Number  on  Dimension¬ 
less  Subcooled  Temperature 

was  increasing  with  increasing  Rcg^i  when  other  parameters  were  fixed.  Therefore, 
based  on  the  Newton’s  law  of  cooling 


q  =  hAT, 

the  subcooled  temperature  difference  (AT)  will  increase  with  increasing  Rtg^i.  This 
corresponds  to  the  trend  indicated  in  Eq.  (4.2)  (note  that  AT*  oc  AT~*). 

4.6  Uncertainty 

The  uncertainties  of  two  phase  pressure  drop  and  Nusselt  number  have  been 
shown  in  Figure  4.3  and  Figure  4.4,  respectively.  The  detailed  uncertainty  anal¬ 
ysis  procedure  is  given  in  Appendix  A.  The  following  are  the  ranges  of  relative 
uncertainties  for  Nu,  X2,  Rtg^i^  and  AP\ 


40 


1.  Uncertainty  of  Nu  =  ±4.6  ~  14.1% 

2.  Uncertainty  of  X2  =  ±0.4  ~  14.9% 

3.  Uncertainty  of  Rcg^i  =  ±8.3  ~  9.6% 

4.  Uncertainty  of  AP  =  ±1.2  ~  180% 

4.7  Discussion 

It  has  been  seen  that  Nusselt  numbers  for  the  chosen  correlations  from  macro¬ 
tubes  have  approximately  3  to  50  times  larger  values  than  the  current  experimen¬ 
tally  determined  Nusselt  numbers.  These  chosen  correlations  represent  a  wide  range 
of  film  condensation  from  stratified  flow  to  annular  flow.  As  mentioned  earlier  in 
this  chapter,  from  the  map  of  flow  pattern  and  the  values  of  the  experimental 
Nusselt  numbers,  the  current  experimental  results  should  be  in  the  spray  flow 
regime  in  which  Nusselt  numbers  should  be  even  greater  than  those  of  annular  flow. 
However,  according  to  Chato’s  recommendation,  the  current  experimental  results 
would  appear  to  be  in  the  stratified  flow  regime.  The  experimental  results  do  not 
agree  with  both  flow  types,  in  particular  the  spray  flow  results.  The  experimental 
Nusselt  numbers  are  much  smaller  than  the  expected  Nusselt  numbers  for  spray 
flow  (approximately  50  times  smaller).  Although  Chato’s  correlation  is  found  to  be 
the  best  of  all  the  correlations  at  predicting  experimental  data,  the  experimental 
data  is  still  a  factor  of  3  less  than  the  expected  values.  In  addition,  it  is  also  found 
that  the  experimental  Nusselt  number  does  not  show  any  dependence  on  Reynolds 
number  (for  the  range  of  Re  values  considered)  or  vapor  velocities,  and  only  depends 
on  the  changes  of  vapor  quality. 

It  has  been  mentioned  previously  that  the  narrow  range  of  Reynolds  numbers 
obtained  in  the  experiment  may  be  the  main  reason  that  an  obvious  trend  of  the 
Nu-Re  relationship  was  not  observed.  As  for  the  dependence  of  Nu  on  the  changes 
of  vapor  quality,  the  increase  in  Nu  with  the  increase  in  Ax  is  due  to  the  increased 
heat  transfer  associated  with  increased  condensation.  Thus,  the  remaining  question 
to  be  answered  is  why  the  experimental  Nusselt  numbers  are  so  small. 
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When  estimating  sizes  of  the  droplets  formed  due  to  condensation,  it  is  appar¬ 
ent  that  droplet  sizes  could  be  on  the  same  order  as  the  inside  diameter  of  the 
microtube.  Bonabian  [42],  Kataoka  et  al.  [43],  Mugele  and  Evans  [44],  and  Wicks 
and  Dukler  [45]  studied  droplet  size  distribution  in  large-tube  annular  two-phase 
flow  without  phase  change  where  droplet  formation  is  mainly  due  to  the  droplet 
entrainment  mechanism.  Kataoka  et  al.  [43]  obtained  the  following  equation  to 
determine  the  average  linear  diameter  of  droplets  Diq: 


D,„  =  0.0031-^fieyq^)  fe)  ,  (4.3 

A>»“|  ^  \PtJ  \PjI 

where  cr  is  the  surface  tension,  and  Rcg  is  the  vapor  Reynolds  number  defined  as 


Re,  =  (4.4) 

Although  the  mechanism  of  droplet  formation  in  the  current  experiment  is 
unclear,  Eq.  (4.3)  may  be  applied  to  determine  the  order  of  magnitude  of  droplet 
sizes.  Using  Eq.  (4.3),  the  droplet  diameter  distribution  for  the  current  experiment 
is  calculated  and  then  plotted  in  Figure  4.6,  where  the  dash  line  represents  the  inside 
diameter  of  the  microtube.  In  this  calculation,  quantities  in  Eq.  (4.3)  including  the 
vapor  velocity  and  all  the  fluid  properties  are  based  on  average  conditions  between 
the  inlet  and  outlet  of  the  microtube.  It  is  seen  from  Figure  4.6  that  droplet 
diameters  are  at  least  on  the  same  order  as  the  inside  diameter  of  the  microtube. 
Thus,  it  seems  reasonable  that  droplet  size  should  be  considered  an  important  factor 
that  affects  the  current  experimental  results. 

Due  to  the  large  droplet  sizes  relative  to  the  inside  diameter  of  the  microtube, 
the  flow  pattern  in  the  current  experiment  should  not  be  categorized  as  spray  flow, 
annular  flow,  or  stratified  flow,  all  of  which  are  based  on  macroscale  concepts.  The 
flow  is  more  likely  an  intermittent  flow  in  which  the  large  liquid  component  fills 
the  entire  cross  section  of  the  microtube  as  shown  in  Figure  2.1.  In  the  current 
case,  the  liquid  component  may  be  considered  to  be  the  result  of  droplet  formation 
along  the  flow  direction  in  the  microtube.  Figure  4.7  schematically  illustrates  an 
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Figure  4.6.  Droplet  Size  Distribution 


expected  flow  pattern.  With  the  alternating  liquid  and  vapor  regions,  the  surface 
tension  force  will  become  important  in  the  current  case. 

The  estimated  shape  of  the  liquid  component  in  the  microtube  may  look  like 
that  shown  in  Figure  4.7.  The  liquid  region  has  a  back  concave  surface  and  a  front 
flat  surface  almost  perpendicular  to  the  flow  direction.  This  is  due  to  the  result  of 
forces  acting  in  the  liquid  including  the  downstream  pressure,  upstream  pressure, 
the  frictional  force,  and  the  surface  tension  force.  The  contact  angles  due  to  the 
surface  tension  for  the  front  surface  and  back  surface  are  assumed  to  be  90°  and 
6  (where  0  ^  90°),  respectively.  This  means  that  the  surface  tension  force  on  the 
back  surface  would  act  solely  against  the  pressure  difference. 

It  was  observed  in  the  experiment  that  the  condensate  liquid  erupted  periodi¬ 
cally  from  the  end  of  the  microtube  into  the  atmosphere.  This  seems  to  support  the 
flow  pattern  theory  proposed  above.  A  force  balance  in  the  streamwise  direction 
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Figure  4.7.  Condensation  Flow  Pattern  Inside  a  Microtube 


on  each  drop  of  liquid,  neglecting  frictional  forces,  gives 

=  ndiCT  cos  0.  (4.5) 

Thus,  the  minimum  pressure  difference  needed  to  overcome  the  surface  tension  force 
is 


AP„  =  —cosO.  (4.6) 

Evaluating  the  surface  tension  of  saturated  water  at  lOO®^  and  1  atm  and  assuming 
a  constant  angle  9  —  45*’,  the  pressure  difference  will  be 


AP,= 


4  X  58.9  X  10  ^A^/mcos45° 
250  X  10~®m 


0.66fcPa. 


(4.7) 


That  is,  at  least  0.5k Pa  pressure  difference  is  needed  to  move  each  drop  of  con¬ 
densate  liquid  in  the  microtube.  Additional  pressure  difference  will  be  needed  to 
overcome  both  frictional  pressure  drop  and  the  momentum  change  pressure  drop 
as  described  in  Chapter  2. 

The  “spitting”  of  fluid  from  the  microtube  outlet  may  be  explained  by  the 
alternating  liquid  and  vapor  regions  within  the  microtube.  In  addition,  the  approx¬ 
imate  IkPa  pressure  difference  required  to  overcome  the  surface  tension  force  and 
to  push  a  liquid  droplet  of  condensate  liquid  from  the  microtube  may  be  unavailable 
at  some  times.  The  local  pressure  on  the  upstream  side  of  the  liquid  region  near 
the  exit  will  have  been  reduced  in  overcoming  the  foregoing  pressure  requirements 
including  those  for  frictional  pressure  drop,  momentum  change  pressure  drop,  and 
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surface  tension  pressure  drop.  Thus,  the  liquid  will  cease  to  flow  temporarily  at  the 
exit  until  the  vapor  pressure  on  the  upstream  side  of  the  liquid  region  is  increased 
by  compression  from  the  adjacent  liquid  region.  When  the  last  liquid  region  is 
ejected  from  the  outlet,  it  will  be  replaced  by  the  adjacent  liquid  region  and  the 
process  will  repeat.  The  result  of  this  flow  pattern  and  process  at  the  exit  is  the 
intermittent  “spitting”  of  liquid  that  was  observed  experimentally. 

The  condensation  heat  transfer  coefficient  with  this  kind  of  flow  pattern  will  be 
smaller  than  that  for  spray  flow,  annular  flow,  or  stratified  flow  since  the  contact 
area  between  the  vapor  and  the  cold  tube  wall  will  be  greatly  reduced  by  the  large 
condensate  drops.  This  is  the  most  likely  explanation  for  the  small  experimental 
Nusselt  numbers. 

It  is  apparent  that  the  two-phase  pressure  drop  inside  microtubes  should  take 
into  account  the  surface  tension  effect.  As  seen  from  Eq.  (4.6),  pressure  drop  due 
to  surface  tension  is  inversely  proportional  to  the  inside  diameter  of  a  microtube. 
Furthermore,  it  is  noticed  that  Eq.  (4.6)  represents  the  pressure  difference  that  is 
used  to  overcome  only  one  drop  of  condensate  liquid.  For  more  than  one  drop  of 
condensate  liquid  within  the  microtube,  the  total  two-phase  condensation  pressure 
drop  should  be  revised  to  be 

AP  =  APj  +  APm  +  NAP,,  (4.8) 

where  N  is  the  number  of  drops  of  condensate  liquid.  Note  that  the  last  term 
in  Eq.  (4.8)  is  only  applicable  in  microtubes  with  liquid  regions  that  completely 
fill  the  cross  section.  In  the  current  experiment,  the  total  pressure  drop  ranged 
from  near  zero  to  65kPa  with  uncertainties  from  ±1.2  ~  180%.  Given  the  level  of 
experimental  uncertainty  and  the  lack  of  data  regarding  the  number  of  condensate 
droplets  in  the  tube,  the  actual  experimental  surface  tension  effect  can  not  be 
ascertained. 
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4.8  Other  Factors  that  may  Affect 
Experimental  Results 

Two  additional  factors  may  have  influenced  the  current  experimental  results. 
These  effects  are  the  result  of  the  presence  of  noncondensible  air  and  high  velocity 
heating. 

In  the  experiment,  tap  water  was  used  as  the  condensing  medium.  Although  it 
is  unclear  how  much  air  is  dissolved  in  the  tap  water,  the  noncondensible  air  will 
certainly  influence  the  condensation  heat  transfer. 

In  a  mixture  of  saturated  steam  and  air,  the  partial  pressure  of  the  steam 
will  be  lower  than  the  total  pressure.  The  reduction  in  pressure  of  the  steam 
will  increase  with  increasing  air  content.  Due  to  this  reduction  of  pressure,  the 
corresponding  saturated  temperature  of  the  steam  will  be  lowered,  resulting  in 
a  smaller  temperature  difference  between  the  steam  and  the  cold  wall.  This  is 
the  reason  that  condensation  heat  transfer  is  degraded  when  noncondensible  air  is 
present.  Stephan  and  Laesecke  [46],  and  Sparrow  et  al.  [47]  investigated  the  forced 
convection  condensation  of  steam  from  a  mixture  with  air.  Their  studies  showed 
that  an  approximately  10%,  20%,  and  30%  reduction  in  condensing  heat  fluxes 
occurred  when  the  partial  pressure  of  air  accounted  for  2%,  5%,  and  10%  of  the 
total  pressure  of  the  mixture,  respectively. 

The  current  experimental  Nusselt  numbers  are  found  to  be  3  to  50  times  lower 
than  those  obtained  from  macroscale  results.  Thus,  noncondensible  air  should 
not  be  a  critical  factor  in  explaining  the  dramatical  decrease  in  condensation  heat 
transfer  in  the  current  experiment. 

Referring  to  Table  4.1,  the  velocity  of  the  inlet  saturated  vapor  is  up  to  270m/s 
in  the  current  experiment.  This  high  velocity  may  be  one  of  factors  that  decreases 
the  condensation  heat  transfer  in  the  microtube  compared  to  macroscale  tubes. 

In  single-phase  flow,  the  heating  effect  causes  the  addition  of  enthalpy  to  the 
fluid  due  to  the  conversion  of  high  velocity  kinetic  energy  of  the  fluid  when  the 
fluid  is  experiencing  a  deceleration  process.  The  effect  associated  with  this  process 
includes  the  addition  of  heat,  and  an  increase  in  temperature  and  pressure  in  the 
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fluid.  The  properties  of  the  fluid  will  also  change  accordingly.  The  Nusselt  number 
in  single  phase  flow  is  related  to  the  heat  transfer  and  temperature  difference  as 


Nu  oc 


Q 

Twi  —  Tf 


Because  both  the  amount  of  heat  (Q)  and  the  temperature  of  the  fluid  (T/)  in  a  high 
velocity  process  include  the  contribution  of  the  high  velocity  heating,  respectively, 
the  Nusselt  number  which  excludes  the  effect  of  the  high  velocity  heating  will 
essentially  be  smaller  than  that  including  the  contribution. 

However,  in  the  condensation  flow  in  a  tube,  the  degree  of  reduction  of  the 
condensation  heat  transfer  coefficient  caused  by  the  effect  of  the  high  velocity 
heating  is  difficult  to  be  estimated.  This  is  primarily  due  to  the  change  of  the 
vapor  velocity  along  the  flow.  For  instance,  the  vapor  velocity  will  be  zero  when  a 
complete  condensation  is  reached,  although  there  has  been  a  high  velocity  of  vapor 
at  the  inlet  of  the  tube. 


CHAPTER  5 


CONCLUSIONS  AND 
RECOMMENDATIONS 

This  thesis  reports  preliminary  experimental  results  of  forced  convection  conden¬ 
sation  of  steam  inside  a  horizontal  microtube.  Conclusions  and  recommendations 
are  given  below. 


5.1  Conclusions 

The  inside  diameter  of  the  tested  microtube  is  the  smallest  so  far  seen  in 
similar  experiments.  It  was  found  that  the  average  experimental  Nusselt  numbers 
of  condensation  (1)  significantly  deviate  from  those  acquired  in  macrotubes,  (2)  do 
not  depend  on  Reynolds  numbers  of  inlet  saturated  vapor  or  inlet  saturated  vapor 
velocities,  and  (3)  only  show  dependence  on  the  changes  of  vapor  quality  across 
the  microtube.  The  significant  deviations  observed  from  the  previously  reported 
macroscale  results  may  be  due  to  a  different  condensation  mechanism  caused  by 
relatively  large  deoplet  sizes  compared  to  the  inside  diameter  of  the  microtube 
tested.  The  two-phase  pressure  should  include  the  contribution  of  surface  tension. 
On  the  other  hand,  it  is  also  noticed  that  the  small  number  of  data  and  other 
limitations  in  the  current  experiment  have  prevented  further  conclusions  to  be 
drawn.  A  large  of  number  of  future  studies  are  needed. 

5.2  Recommendations 

In  view  of  limitations  involved  in  the  current  experiment,  some  suggestions  for 
further  studies  are  offered  as  follows. 

1.  Include  a  wider  range  of  sizes  of  microtubes  in  order  to  investigate  effects  of 
geometries  on  condensation.  This  may  include  finding  the  sizes  at  which  con- 
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densation  Nusselt  numbers  or  pressure  drops  start  to  deviate  from  macroscale 
results. 

2.  Test  different  kinds  of  fluid  to  study  effects  of  fluid  properties  on  condensation. 

3.  Increase  the  range  of  Reynolds  numbers  to  better  observe  the  expected  Re 
effects. 

4.  Study  effects  of  surface  roughness  on  condensation. 

5.  Investigate  effects  of  orientation  of  microtubes  because  correlations  for  annu¬ 
lar  flow  condensation  available  in  the  literature  are  also  applicable  to  non¬ 
horizontal  macrotubes. 

6.  Improve  experimental  accuracy.  This  can  be  realized  by  reducing  number  of 
intermediate  steps  that  are  used  to  determine  the  average  Nusselt  numbers. 
For  instance,  using  copper  instead  of  stainless  steel  in  the  current  experiment 
is  expected  to  improve  the  accuracy  of  experimental  average  Nusselt  numbers 
because  the  effects  of  thermal  conduction  in  the  microtube  walls  may  be 
ignored  due  to  the  much  higher  thermal  conductivity  of  copper  than  that 
of  stainless  steel.  Additionally,  the  accuracy  can  be  increased  by  directly 
measuring  some  quantities  instead  of  making  assumptions  and  simplifications 
for  them.  For  instance,  directly  measuring  the  inlet  vapor  quality  would 
improve  the  experimental  accuracy. 

7.  Visualize  the  process  of  condensation  to  obtain  a  better  understanding  of 
the  formation  of  condensate  and  the  distribution  of  condensate  or  vapor  in 
microtubes.  To  do  so,  using  microchannels  instead  of  microtubes  may  be  a 
better  scheme. 


APPENDIX  A 


UNCERTAINTY  ANALYSIS 

For  a  result,  F,  which  is  a  function  of  a  set  of  variables  (xi,a;2,-  •  the 

absolute  uncertainty  of  F  is 


This  is  the  basic  working  equation  for  an  uncertainty  analysis.  In  this  equation, 
is  called  the  sensitivity  coefficient  which  can  indicate  how  sensitive  the  variation 
of  the  function,  F ,  to  the  variation  of  the  variable  x,-  will  be,  represents  the 

contribution  of  the  variable  x,-  to  the  total  uncertainty  (dF)  of  the  function  F,  and 
Sxi  is  the  uncertainty  interval.  The  uncertainty  interval,  5x,-,  may  usually  be  set  to 
be  half  the  smallest  scale  division  of  the  instrument.  Relative  uncertainty  of  F  is 
defined  as 

^  X  100%. 
t 

Uncertainties  of  current  experimental  results  are  obtained  using  an  uncertainty 
reduction  computer  program  that  is  made  based  on  the  Moffat  method  [48].  The 
key  point  in  this  program  is  to  find  the  sensitivity  coefficients  by  using  the  central- 
difference  definition  of  a  partial  derivative,  that  is, 

/^(xi,  X2,  •  •  -,  X,-  -f  £,  •  •  •,  Xn)  —  F"(xi,  X2,  •  •  -,  X,  —  £,  •  •  •,  X^) 

dxi  '-»•<»  2e  ’ 

where  the  incremental  variable,  e,  may  be  set  to  a  sufficiently  small  number  to 
obtain  accurate  enough 

In  the  uncertainty  reduction  program,  e  =  10~®Xi  for  all  the  15  variables  which 
are  used  to  determine  Nusselt  number  (Nu),  vapor  quality  (X2)  at  the  outlet  of 
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the  microtube,  Reynolds  number  {Rcg^i)  of  inlet  saturated  vapor,  and  two-phase 
pressure  drop  (AP).  The  determination  of  measurement  uncertainty  intervals,  ^x,-, 
of  the  15  variables  are  based  on  data  provided  by  manufacturers  of  the  instruments. 
These  uncertainty  intervals  are  listed  in  Table  A.l. 

The  following  is  a  typical  set  of  uncertainty  data  calculated  by  using  the  above 
uncertainty  reduction  program: 

Raw  experimental  data  are 

1.  di  =  0.25mm 

2.  do  =  0.79mm 

3.  L  =  47.90mm 

4.  D  =  4.93mm 

,  5.  =  23.3°C 

6.  T^,o  =  52.5°^ 

7.  Zpht  =  123.2°C 

8.  Tw,o,i  =  86.3°C 

9.  T^i,o,2  =  90.3°C 


Table  A.l.  Uncertainty  Intervals 


No. 

Xi 

Sxi 

1  No. 

Xi 

Sxi 

m 

Sdi 

0.02(mm) 

mm 

o.i%r^i,o,2  +  o.6(°c) 

m 

Sdo 

0.02(mm) 

10 

IHIlii 

o.i%t^„/,o,3  +  o.eC’C') 

IBl 

5L 

0.02(mm) 

11 

0.1%T^i,o,4  +  0.6{°C) 

4 

5D 

b.02(mm) 

12 

SPo 

0.ee7{kPa) 

0.i%r,,,;  -t-  0.6(°C') 

13 

2.~59{kPa) 

im 

o.i%r,„,,-i-o.6("C') 

14 

SVod 

2%Vcd{ml !  min) 

^Tgpht 

0.1%Tspht  +  0.0{°C) 

15 

sv^ 

2%Vcd{ml !  min) 

8 

0.1%Tw.o,i  +  0.6("C) 
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10.  T^i,o,3  =  %.2°C 

11.  T^i,oa  = 

12.  Po  =  M.mkPa 

13.  Papht  =  32.19A:Pa 

14.  Vcd  =  0.3400m//mm 

15.  14;  =  1.7800m//mm 

The  corresponding  uncertainties  for  a;2,  TVu,  and  AP  are 

1.  =  2395.6  ±  210.7 

2.  =  0.7328  ±0.0114 

3.  =  3.7670  ±  0.3061 


4.  AP  =  6.7592  ±  9.0678itPa 


APPENDIX  B 


EXPERIMENTAL  DATA 


The  following  table  includes  46  sets  of  experimental  data: 


No. 

i 

2 

3 

4 

5 

6  7 

8 

22.8 

22.8 

22.8 

iBn 

23.2 

T^A°C) 

48.7 

50.0 

mm 

\mm 

51.3 

Esn 

T^i,oA°C) 

85.5 

84.2 

84.2 

84.8 

82.6 

84.1 

84.8 

mmm 

T^i,oA°C) 

mm 

89.1 

mm 

89.1 

89.4 

87.3 

Tn,l,oA°C) 

88.9 

88.2 

\mm 

88.0 

89.4 

88.0 

86.5 

Tu>1,o,4{^C) 

91.8 

bIKM 

91.1 

KIMI 

92.2 

92.0 

91.8 

Zphti°C) 

118.3 

urn 

115.3 

Vcd(Tnl  /  min) 

0.098 

0.116 

0.166 

lillrAl 

itiwaa 

QQI 

|[QQ| 

Po{kPa) 

maam 

84.70 

84.72 

84.73 

Pspht{kPa) 

11.57 

11.16 

11.16 

16.62 

15.96 

No. 

9 

10 

11 

m 

13 

14 

16 

T^A^C) 

22.4 

mM 

22.5 

22.9 

22.9 

22.9 

23.3 

mmm 

T^A°C) 

mm 

mm 

46.7 

49.4 

53.7 

46.5 

52.5 

mm 

T^i,oA^C) 

84.7 

80.9 

88.0 

83.7 

79.8 

Tv.i,oA°C) 

m 

89.0 

86.3 

88.8 

86.2 

90.8 

88.2 

%uoA°C) 

KIIm 

89.0 

87.3 

mam\ 

89.4 

86.7 

90.2 

88.1 

T^i,oa{°C) 

90.0 

88.5 

mm 

90.3 

89.2 

86.1 

89.5 

88.0 

Tsphii^'C) 

130.7 

ma 

116.0 

118.2 

116.7 

Vcdiml/min) 

0.234 

0.234 

0.202 

0.340 

0.356 

Ki,(m//min) 

1.120 

3.630 

UMil 

2.390 

3.6300 

1.780 

3.790 

PAkPa) 

84.66 

mm 

84.70 

1^ 

84.72 

84.71 

84.75 

P sphtikP 

75.14 

71.08 

61.71 

15.21 

mm 

md 

53 


No. 

17 

18 

19 

20 

21 

22 

23 

T^A°C) 

23.1 

Mm 

Hi 

22.9 

KfeliM 

23.1 

T^A°C) 

33.3 

45.2 

IBiCT 

i^a 

53.4 

mmm 

Twl,o,\{^C) 

75.5 

89.6 

87.6 

85.7 

89.2 

mm 

f^i,oa{°C) 

86.8 

91.0 

90.1 

90.6 

88.1 

ireafeM 

T^i,oA°C) 

86.8 

91.0 

90.1 

90.6 

88.1 

89.2 

TwI,o,a{^C) 

87.2 

88.6 

85.8 

89.3 

89.2 

89.0 

88.4 

89.1 

Tspht{°C) 

116.7 

110.7 

110.9 

117.1 

125.7 

120.3 

120.1 

120.3 

Vcd{ml/min) 

0.356 

0.056 

0.128 

0.283 

0.193 

0.202 

0.270 

Biifcfii 

VA'fnl'  f  min) 

lUjUjlll 

imjj 

PoikPa) 

HI 

84.73 

84.69 

84.71 

84.71 

84.73 

84.70 

Pspht{kPa) 

24.81 

28.20 

No. 

28 

TU^C) 

23.1 

22.3 

22.5 

EOH 

T^A°C) 

"48.7 

mm 

49.1 

54.6 

53.6 

mm 

47.2 

88.2 

89.6 

87.9 

90.0 

87.0 

90.5 

82.5 

TwI,o,2\C) 

87.8 

88.5 

91.5 

88.0 

85.3 

88.30 

iaoi 

90.5 

92.5 

90.7 

90.6 

91.2 

88.9 

T^I,o,4{‘'C) 

88.4 

liBII 

92.3 

MM 

Zpkti^C) 

112.1 

l^ll 

■iiiia 

■KWIM 

VRISIEI 

Vcdiml/min) 

0.158 

DBEII 

mm 

mm\ 

BUM 

BliW'J 

Vwiml/min) 

0.299 

0.380 

0.947 

0.792 

1.660 

1.550 

1.660 

3.630 

PoikPa) 

84.72 

84.73 

84.73 

84.75 

84.75 

84.76 

84.75 

84.76 

Pspht{kPa) 

20.49 

72.44 

70.00 

57.20 

45.53 

54.48 

51.06 

53.78 

54 


No. 

33 

34 

35 

36 

37 

38 

39 

mm 

21.9 

21.1 

11^ 

22.4 

22.6 

22.6 

22.6 

BarjM 

T^A^C) 

38.2 

56.5 

wm 

39.6 

53.4 

52.1 

55.0 

EpmE 

TwI,oa{^^) 

78.9 

87.2 

80.0 

79.7 

87.7 

90.8 

83.9 

T^i,oA°C) 

ion 

84.0 

mmm 

89.1 

91.2 

T^i,oA°C) 

11^ 

QISIB 

88.2 

'90.7 

91.3 

Tyjl^o,A{^C) 

87.8 

89.7 

88.5 

87.8 

90.6 

90.3 

89.1 

123.1 

136.9 

123.8 

130.4 

128.1 

Vcd{mllmin) 

0.356 

0.560 

0.213 

0.425 

0.425 

VAiT^ilfnin) 

IQgg 

2.260 

QQjl 

0.721 

2.780 

5.360 

PoikPa) 

84.77 

84.77 

84.77 

84.77 

Pspht{kPa) 

52.49 

98.39 

96.33 

96.14 

43.22 

38.54 

78.95 

71.60 

42 

43 

46 

T^A'C) 

22.8 

22.6 

22.6 

22.8 

T^A°C) 

55.2 

49.5 

mmm 

37.6 

T..i,oA°C) 

86.7 

87.6 

83.3 

79.5 

80.2 

T^i,oA°C) 

88.1 

89.1 

84.1 

87.5 

90.7 

90.9 

88.0 

89.4 

90.5 

89.5 

1^ 

88.0 

88.1 

88.6 

lEMSH 

133.7 

133.6 

■ItTilM 

129.4 

0.521 

0.643 

■W«M 

Em 

2.260 

HH 

KHatllM 

7.360 

7.360 

3.630 

PoikPa) 

84.76 

Riwa 

84.78 

84.79 

omi 

92.33 

101.70 

97.30 

95.66 

imi 
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